Method of magnetic resonance arteriography using contrast agents

ABSTRACT

The present invention is a method and apparatus for providing preferential enhancement of an artery of interest relative to adjacent veins and background tissue. The method and apparatus adapts the timing of a maximum or substantially elevated rate of infusion to correlate with the collection of image data corresponding to the center of k-space. The technique and apparatus temporally correlates the timing of a maximum or substantially elevated rate of infusion and the mapping of k-space according to the location of the artery of interest, the size of the artery of interest, the physical condition of the patient, the time delay due to the configuration of the contrast agent delivery system, and/or the type of pulse sequence employed by the imaging apparatus. Adapting the timing of a maximum or substantially elevated rate of infusion to correlate with the collection of image data corresponding to the center of k-space provides a period of a maximum or substantially elevated contrast concentration in the artery of interest relative to adjacent veins during collection of at least a portion of the image data corresponding to the center of k-space.

It is hereby noted that work associated with reducing this invention to practice was supported, in part, by the National Institutes of Health, under Grant No. HL 46384. The United States Government may have certain rights in this invention.

This application is a continuation of application Ser. No. 08/378,384, filed Jan. 25, 1995 now U.S. Pat. No. 5,553,619 which is a continuation-in-part of application Ser. No. 08/071,970, filed on Jun. 7, 1993 now U.S. Pat. No. 5,417,213.

BACKGROUND OF THE INVENTION

This invention relates to a method of, and apparatus for use in, magnetic resonance angiography to detect, diagnose, and treat arterial diseases and injuries. Arterial diseases and injuries are common and often have severe consequences including death. Imaging arteries serves to detect and characterize arterial disease before these consequences occur.

A conventional method of arterial imaging includes inserting a catheter into the artery of interest (the artery under study) and injecting radiograpic contrast, for example, an iodinated contrast, while taking radiographs of the artery. Radiographs are commonly referred to as X-rays. In this technique, the contrast remains in the arteries for a few seconds during which the arteries appear distinct from both the veins and background tissue in the radiographs.

Although a catheter-based contrast arteriography technique generally provides high quality arterial images, there is a risk of arterial injury or damage by the catheter and its insertion. There may be thrombosis, dissection, embolization, perforation or other injury to the artery itself. Furthermore, such a technique may result in a stroke, loss of a limb, infarction or other injury to the tissue supplied by the artery. In addition, hemorrhage at the catheter insertion or perforation sites may require blood transfusions. Moreover, kidney failure and brain injury may result from the toxic effects of the X-ray contrast.

More recent techniques of arterial imaging are based upon detecting the motion of the blood within the arteries and/or veins. These techniques involve employing magnetic resonance imaging (MRI) to image moving blood distinct from stationary background tissues. (See, e.g., Potchen, et al., eds., "Magnetic Resonance Angiography/Concepts and Applications", Mosby, St. Louis, 1993; the text of which is incorporated herein by reference). Such techniques do not necessitate catheter insertion into the artery. These techniques are commonly known as 2D time-of-flight, 3D time-of-flight, MOTSA, magnitude contrast, phase contrast, and spin echo black blood imaging.

With pre-saturation pulses it is possible to primarily image blood flowing in one direction. Since arteries and veins generally flow in opposite directions, these pre-saturation pulses allow preferential visualization of the arteries or the veins. Because these techniques depend upon blood motion, the images are degraded in patients who have arterial diseases which decrease or disturb normal blood flow. Such types of arterial diseases that decrease or disturb normal blood flow include aneurysms, arterial stenoses, arterial occlusions, low cardiac output and others. The resulting lack of normal blood flow is particularly problematic because it is those patients with disturbed blood flow in whom it is most important to acquire good quality arterial images.

A related magnetic resonance imaging technique relies upon differences in the proton relaxation properties between blood and background tissues. (See, e.g., Marchal, et al., in Potchen, et al., eds., supra, pp. 305-322). This technique does not depend upon steady blood in-flow. Instead, this magnetic resonance imaging technique involves directly imaging the arteries after administering a paramagnetic contrast agent. Here, after administering the contrast agent, it is possible to image arteries directly based upon the blood relaxation properties. This technique overcomes many of the flow related problems associated with magnetic resonance imaging techniques which depend upon blood motion.

Several experts have performed magnetic resonance arterial imaging using intravenous injection of gadolinium chelates (paramagnetic contrast agents). These experts have reported their results and conclusions. In short, these results have been disappointing and, as a result, the use of gadolinium for imaging arteries has not been adopted or embraced as a viable arterial imaging technique. The images using this technique are difficult to interpret because the gadolinium tends to enhance both the arteries and the veins. Since the arteries and veins are closely intertwined, it is extremely difficult to adequately evaluate the arteries when the veins are visible. Further, the difficulty in interpretation is exacerbated as a result of contrast leakage into the background tissues.

As a result, there exists a need for an improved method of magnetic resonance angiography which provides an image of the arteries distinct from the veins and which overcomes the limitations of other techniques. Further, there exists a need for an apparatus which facilitates providing an image of the arteries distinct from the veins and which may be implemented in overcoming the limitations of other techniques.

SUMMARY OF THE INVENTION

In one aspect, the present invention is a method of imaging an artery in a patient using magnetic resonance imaging. The method includes administering magnetic resonance contrast agent to the patient prior to and/or during collection of image data. The contrast may be administered by intravenous infusion at a rate of infusion sufficient to provide a substantially elevated or maximum concentration of the contrast agent in the artery during collection of image data representative of a center of k-space.

In a preferred embodiment, the administration of the contrast agent includes temporally correlating a substantially elevated or maximum rate of infusion of the contrast agent with the collection of image data representative of the center of k-space. The rate of infusion and collection of image data are correlated in accordance with a delay time in a delivery system, a location of the artery and/or a physical condition of the patient.

In another preferred embodiment, the method according to this aspect of the invention further includes administering the contrast agent at a substantially elevated or maximum rate of infusion about 10 to about 30 seconds before collection of image data representative of the center of k-space. Where the contrast agent is paramagnetic, a maximum rate of infusion may be greater than 0.0015 Liters/Kg-sec² divided by the relaxivity.

Further, in another preferred embodiment, the method may include temporally correlating a period of a substantially elevated rate of infusion with the collection of image data representative of the center of k-space in accordance with a size of the artery so that the concentration of the contrast agent in the artery is substantially greater than veins adjacent to the artery during collection of image data representative of the center of k-space. In this preferred embodiment, the contrast agent may be administered at an infusion rate and for a period which provides a substantially elevated or maximum concentration of the magnetic resonance contrast agent in the artery during more than 50% of the period of collecting image data representative of the center of k-space. The contrast may also be administered to provide a substantially elevated or maximum concentration of the magnetic resonance contrast agent in the artery for a period of between about 50% to about 85% of the time of collecting image data representative of the center of k-space.

In one preferred embodiment, the step of administering the contrast agent to the patient further includes administering a paramagnetic contrast agent having a relaxivity. The paramagnetic contrast agent may be administered at a rate of infusion sufficient to provide a concentration of the paramagnetic contrast agent in the artery, during collection of image data representative of a center of k-space, of greater than 2.9/sec*relaxivity.

In another preferred embodiment, the techniques of this invention is implemented using a mechanical injector which is spring-loaded, pneumatic powered, or electrically powered. The mechanical injector may include a non-magnetic spring to pressurize the contrast agent for infusion into the patient. Further, the mechanical injector may be adapted to receive a vessel containing a gadolinium chelate.

In another aspect, the present invention is a method for administering a magnetic resonance contrast agent into a vein of a patient to enhance a magnetic resonance image of an artery of the patient. The method includes administering the contrast agent into the vein at an infusion rate and for a period of time sufficient to provide a maximum or substantially elevated concentration of contrast in the artery during collection of at least a portion of the image data which corresponds to a center of k-space.

In a preferred embodiment of this aspect of the invention, the method further includes administering a paramagnetic contrast agent at a maximum rate of infusion which is greater than 0.0015 Liters/Kg-sec² divided by the relaxivity.

In another preferred embodiment, the method may further include temporally correlating a maximum rate of infusion of the contrast agent with the mapping of k-space in accordance with a delay time in a delivery system, a location of the artery or a physical condition of the patient.

The paramagnetic contrast agent may be administered at a rate of infusion sufficient to provide a concentration of the paramagnetic contrast agent in the artery of greater than 2.9/sec*relaxivity during collection of image data representative of a center of k-space.

The method of this aspect of the invention may also be implemented using a mechanical injector which is spring-loaded, pneumatic powered, or electrically powered. The mechanical injector may be adapted to receive a vessel containing a gadolinium chelate. Further, the mechanical injector may include a non-magnetic spring to pressurize the contrast agent for infusion into the patient.

In yet another aspect, the present invention is a magnetic resonance compatible apparatus for administering a magnetic resonance contrast agent into a vein of a patient to enhance a magnetic resonance image of an artery of the patient. The apparatus includes a spring-powered device, having a non-magnetic spring, for administering the contrast agent into the vein at an infusion rate and for a period of time sufficient to provide a substantially elevated concentration of the contrast agent in the artery during collection of at least a portion of the image data which corresponds to a center of k-space. The non-magnetic spring may be manufactured from eljaloy, inconel or other non-magnetic material.

The magnetic resonance contrast agent, in this aspect of the invention, may be contained in a syringe having an orifice at a distal end and a plunger for engaging the contrast agent. The syringe may be coupled to the spring-powered device such that the spring-powered device engages the plunger of the syringe to administer the contrast agent into the vein at an infusion rate within the range of about 4 ml/minute to about 600 ml/minute.

In a preferred embodiment, the apparatus includes a rate adjustment mechanism coupled to the plunger of the syringe to adjust the infusion rate of the contrast agent. Further, the apparatus may also include a fluid flow restrictor coupled to the orifice of the syringe to control, in combination with the non-magnetic spring, the infusion rate of the contrast agent. The fluid flow restrictor may be a needle, a precision orifice, or tubing having a narrow calibre.

The apparatus may also include a flow rate indicator to visually or audibly indicate the infusion rate of the contrast.

In still yet another aspect, the present invention is an apparatus for administering a magnetic resonance contrast agent into a vein of a patient to enhance a magnetic resonance image of an artery of the patient. The apparatus includes an infusion mechanism, including infusion means for matching the infusion rate with the mapping of at least a portion of k-space. The infusion means infuses the contrast agent at a rate and duration sufficient to provide a maximum or elevated concentration of the contrast agent in the artery during collection of at least a portion of the image data which corresponds to a center of k-space.

Finally, in another aspect, the present invention is an apparatus for injecting a magnetic resonance contrast agent into a patient, including an injection mechanism to match a substantially elevated or maximum contrast agent concentration in the artery relative to adjacent veins with a mapping of k-space by controlling the rate of flow of the contrast agent such that a substantially elevated rate of flow correlates with the collection of image data corresponding to a center of k-space.

In a preferred embodiment, the injection mechanism may temporally correlate the substantially elevated or maximum rate of infusion with the collection of image data corresponding to the center of k-space. This correlation may be in accordance with a location of the artery and/or the time delay inherent in the infusing apparatus.

The injection mechanism may provide a substantially elevated or maximum rate of infusion of the contrast agent about 10 to about 40 seconds before beginning collecting image data which corresponds to the center of k-space.

In a preferred embodiment, the injection mechanism temporally correlates a period of a substantially elevated or maximum rate of infusion of the contrast agent with the collection of image data corresponding to the center of k-space in accordance with a size of the artery. This provides for a situation where the concentration of the contrast agent in the artery is substantially greater than adjacent veins while collecting image data which corresponds to the center of k-space. In another preferred embodiment, the injection mechanism includes means for infusing the contrast agent into the patient at a rate and duration of infusion which provides a substantially elevated or maximum concentration of the magnetic resonance contrast agent in the artery for at least 50% of the time of collecting image data corresponding to the center of k-space.

The present invention overcomes the limitations of other techniques by injecting magnetic resonance contrast agents at a sufficient rate, at a selected time relative to the collection of image data, and for an appropriate duration in such a manner that the contrast level in the arteries is higher than that in surrounding veins and background tissue during collection of image data. The injection may be intravenously in a vein remote from the artery of interest. Intravenous injection eliminates the risks associated with arterial catheterization. In the present invention, the high level of arterial contrast permits directly imaging the arterial lumen, analogous to conventional arteriography. Moreover, using a magnetic resonance pulse sequence which is not as sensitive to motion and by relying on image contrast related to differences in T1 relaxation rather than the in-flow effect, a reduction in the flow artifacts associated with phase contrast or magnitude contrast (time-of-flight) magnetic resonance angiography is observed.

In short, the present invention is, in comparison or relative to other techniques, a method of magnetic resonance angiography which combines several of the advantages of catheter-based contrast arteriography with the advantages of magnetic resonance imaging while substantially eliminating the disadvantages of each.

BRIEF DESCRIPTION OF THE DRAWINGS

In the course of the detailed description of preferred embodiments to follow, reference will be made to the attached drawings, in which:

FIG. 1 illustrates longitudinal relaxation time (T1) of blood as a function of injection imaging time and total paramagnetic contrast dose for a compound with a relaxivity of 4.5/millimolar-second;

FIG. 2 illustrates calculated magnetic resonance signal intensity as a function of flip angle for 5 different longitudinal relaxation times (T1) assuming a spoiled, 3D volume acquisition with TR equal to 25 msec and TE<<T2*;

FIGS. 3, 4, 5A and 5B, and 6A-C are block diagram representations of mechanical infusion devices and configurations, according to the present invention;

FIG. 7 is a block diagram representation of a manual injection configuration, according to the present invention;

FIGS. 8A-C illustrate typical coronal maximum intensity projection (MIP) collapse images obtained (FIG. 8A) prior to injection of gadopentetate dimeglumine, (FIG. 8B) dynamically during intravenous injection of gadopentetate dimeglumine, 0.2 millimoles/kilogram over 5 minutes, and (FIG. 8C) immediately following injection of gadopentetate dimeglumine;

FIG. 9 illustrates region of interest analysis averaged for 3 patients who had pre-infusion, dynamic infusion, immediate post infusion and delayed 3D FT imaging. This figure shows that there is a short window, during contrast infusion, when the aorta signal intensity (solid squares) is higher than that of the IVC (open squares) and background tissues, fat (diamonds) and muscle (triangles);

FIG. 10A is an illustrative example of a magnetic resonance image of a patient with an abdominal aortic aneurysm. The magnetic resonance angiography ("MRA") depicts the aneurysmal aorta and aneurysmal common iliac arteries as well as severe stenoses of the right external iliac (curved arrow) and inferior mesenteric (straight arrow) arteries and a mild stenosis of the left common iliac artery; and

FIG. 10B illustrates a digital subtraction angiogram of the aneurysmal aorta and aneurysmal common iliac arteries as well as severe stenoses of the right external iliac (curved arrow) and inferior mesenteric (straight arrow) arteries and a mild stenosis of the left common iliac artery of FIG. 10A.

DESCRIPTION OF A PREFERRED EMBODIMENTS

Briefly, the present invention is a technique and apparatus for providing preferential enhancement of an artery of interest relative to adjacent veins and background tissue by adapting the timing of a maximum or substantially elevated rate of infusion to correlate with the collection of image data representative of the center of k-space. The technique and apparatus temporally correlates the timing of a maximum or substantially elevated rate of infusion and the mapping of k-space according to the location of the artery of interest, the size of the artery of interest, the physical condition of the patient, the time delay due to the configuration of the contrast agent delivery system, and/or the type of pulse sequence employed by the imaging apparatus. Adapting the timing of a maximum or substantially elevated rate of infusion to correlate with the collection of image data corresponding to the center of k-space provides a period of a maximum or substantially elevated contrast concentration in the artery of interest relative to adjacent veins during collection of at least a portion of the image data representative of the center of k-space.

In particular, in one aspect, the present invention is a method of magnetic resonance angiography wherein the magnetic resonance imaging data is collected over a period with simultaneous controlled intravenous injection of a magnetic resonance contrast agent. The magnetic resonance contrast agent is preferably injected into a patient, for example, a human or other animal, substantially throughout the period of imaging in a controlled manner, i.e., injected at a controlled rate over the period of imaging. In a preferred embodiment, the magnetic resonance contrast agent is administered as a steady and continuous infusion in a vein which is remote from the artery of interest (i.e., the artery under study).

The method and apparatus of this invention is especially useful in imaging peripheral arteries. A peripheral artery is an artery that is outside the brain. Since the blood-brain barrier prevents redistribution of the contrast agent from the capillary bed in the brain, this invention is more effective in imaging arteries found in a patient's body which are below the neck.

Magnetic resonance contrast agents employed in the present invention are well known in the art, and are disclosed in, for example, U.S. Pat. Nos. 5,141,740; 5,078,986; 5,055,288; 5,010,191; 4,826,673; 4,822,594; and 4,770,183, which are incorporated herein by reference. Such magnetic resonance contrast agents include many different paramagnetic contrast agents, for example, gadolinium compounds. Gadopentetate dimeglumine and gadoteridol are paramagnetic gadolinium chelates that are readily available, and which rapidly redistribute into the extracellular fluid compartment. Other gadolinium compounds are acceptable, and may have a higher relaxivity, more rapid redistribution into the extracellular fluid compartment, and greater and faster extraction in the capillary bed. It should be noted that, contrast agents that are extracted or degrade in or near the capillary bed are preferred for the present invention.

In one preferred embodiment, the injected magnetic resonance contrast agent should be sufficiently small to rapidly redistribute into the extracellular fluid compartment in the systemic capillary bed, or the contrast agent should be actively extracted from the circulation in the capillary bed distal to the artery of interest, or both. Under these circumstances, the artery (or arteries) of interest contains a high concentration of contrast and the vein (or veins) adjacent to the artery (or arteries) of interest possesses a low contrast concentration. Further, under these circumstances, the relationship of artery-to-venous contrast concentration is substantially maintained over the period of contrast injection.

By matching the duration of the injection with the time required for a longitudinal relaxation time (T1) weighted magnetic resonance image data set, it is possible to view the arteries distinct from the veins. Further, by injecting the contrast at a sufficient rate, the longitudinal relaxation time of the arterial blood may be made sufficiently short when compared to that of the background tissues. As a result, the image of the arteries is distinct from background tissue as well.

As mentioned above, the magnetic resonance contrast agent is administered to the patient, for example, a human or other animal, via intravenous infusion, i.e., injection into a patient at a controlled rate over a period of time. In a preferred embodiment, the period of infusion of magnetic resonance contrast agent may be a substantial portion of the time during which image data is being collected for a magnetic resonance image. A substantial portion of the data collection time is a majority of the time and should include the period of time during which the center of k-space is acquired. A substantial portion of the data collection may be at least 75% of the collection period, and preferably greater than 85% of the collection period.

In another preferred embodiment, a maximum or elevated rate of infusion is correlated with collecting of image data representing the center of k-space. That is, a maximum or elevated rate of injection of contrast occurs while mapping the center of k-space. The center of k-space refers to the low spatial frequency MR image data. In a preferred embodiment, the timing of the infusion of the magnetic resonance contrast agent (e.g., gadolinium), however, is compensated to account for the time which may be necessary for the contrast agent to circulate from the site of injection, through the heart and lungs, and into the artery of interest. Under this circumstance, in order to more accurately correlate a maximum or elevated rate of infusion of the contrast agent with collection of image data representative of the center of k-space, it may be necessary for the period of a maximum or substantially elevated rate of infusion of contrast to be forwardly shifted in time relative to the acquisition or collection of image data representing the center of k-space.

Correlating the maximum or substantially elevated rate of infusion with the mapping of k-space accounts for the time delay due to the contrast agent delivery system (e.g., the length of catheter which delivers the contrast agent) and/or the time-delay due to the time required for the gadolinium to circulate from the site of injection, through the body, and into the artery of interest. Further, correlating a maximum or elevated rate of infusion with the mapping of the center of k-space provides for a period of maximum or elevated arterial contrast in the artery of interest to occur during collection of image data corresponding to the center of k-space (which in most pulse sequences is collected during the middle of a scan).

Typically, there is a delay of about 10-40 seconds from the time of administering the contrast agent at a maximum or substantially elevated rate to the realization of a maximum or elevated arterial contrast in the artery of interest relative to adjacent veins. As a result, the amount of time necessary to align the collection of data representing the center of k-space with a period of maximum or elevated arterial contrast is about 10 to 40 seconds after initiating maximum or substantially elevated rate or infusion. That is, the maximum or substantially elevated rate of infusion should occur about 10 to 40 seconds prior to acquiring image data representative of the center of center of k-space to "compensate" for the contrast delivery system and the circulation time of the contrast agent. In this regard, the maximum rate of infusion of the contrast agent is correlated with the collection of image data which corresponds to the center of k-space in that a period of maximum or elevated arterial contrast occurs during the mapping of the center of k-space.

In some instances, the pulse sequence of the imaging apparatus may be designed such that data representing the center of k-space is collected at the beginning of a scan. Under this circumstance, the above described technique may be implemented by injecting contrast agent prior to the image scanning process (collection of image data) so that the injection of the contrast agent is properly correlated to the acquisition of data corresponding to the center of k-space and a period of maximum or elevated arterial contrast coincides with that data acquisition.

In one preferred embodiment, the timing of the infusion of the contrast agent may be further adjusted according to the physical condition of the patient. That is, the timing of the administration of the contrast agent to, for example, correlate a maximum or elevated rate of infusion with the mapping of k-space, may be adjusted according to a patient's physical condition which may affect the amount of time for contrast to circulate from the site of infusion (typically an arm vein) to the artery of interest. Although this time interval will typically be about 10-40 seconds, patients with cardiac disease or under sedation may have a slower flow, which may suggest a time interval of about 30 seconds or greater.

In contrast, young healthy patients, and especially patients who have recently exercised, will have a faster flow, which may suggest a time interval of less than 20 seconds. Faster flow is also seen in patients with arterio-venous fistulas. The circulation time will also be decreased in patients who have central lines that allow injection of contrast close to or directly into the heart.

In another preferred embodiment, the timing of a maximum or substantially elevated infusion rate may be further selected or controlled, based on the location of the artery of interest relative to the injection site and/or the patient's heart, to provide a period of maximum or elevated arterial contrast concentration in the artery of interest relative to adjacent veins during collection of image data which corresponds to the center of k-space. In this regard, the amount of time between a maximum or elevated rate of injection and a maximum or elevated arterial contrast is realized in the artery of interest may be longer if the artery of interest is far from the heart, for example, in the foot, or if the site of infusion is a peripheral vein such as in the hand. Thus, in this preferred embodiment, the correlation of the maximum rate of injection with the collection of image data representative of the center of k-space includes adapting the time between the maximum rate of injection and the mapping of the center of k-space to reflect the relative location of the artery of interest.

In another preferred embodiment, the correlation between a maximum or elevated rate of injection and the collection of image data representative of the center of k-space may be adjusted or controlled in accordance with the size of the artery of interest. This is because the varying image contrast during image acquisition creates artifacts. If the arteries have a high concentration of contrast for only a short period of time when the center of k-space is acquired, then the artifacts will be larger, making it possible to evaluate only the largest of arteries. In order to see small arteries free of artifacts, it is necessary to have a high concentration of gadolinium in the artery for a substantial fraction of the image acquisition time, including the period corresponding to acquisition of the center of k-space. In this regard, the duration for which the magnetic resonance contrast agent is administered at a maximum or elevated rate during collection of image data corresponding to the center of k-space may be adapted according to the size of the artery of interest. This provides for a period of overlap between a maximum or elevated concentration of contrast in the artery of interest and mapping of the center of k-space (i.e., overlap correlation factor). The overlap correlation factor may be defined as the amount of time during which an elevated or maximum arterial contrast exists in the artery of interest relative to the collection of image data representing the center of k-space. The overlap correlation factor may also be defined as the duration of overlap between the collection of image data representative of the center of k-space and a maximum or elevated rate of injection of the contrast agent. Where the artery of interest is relatively large (e.g., the aorta), an elevated or maximum injection rate of the contrast agent may overlap with the collection of image data representative of the center of k-space for a shorter period of time. In contrast, where the artery is relatively small (e.g., a renal artery), the contrast agent may be administered at an elevated or maximum rate for a longer period of time during the mapping of the center of k-space.

Further, when imaging larger arteries, in a preferred embodiment, an elevated or maximum concentration of contrast agent in the artery of interest is provided for at least 20% of the time during which image data corresponding to the center of k-space is collected; and preferably, an elevated or maximum concentration of contrast in the artery of interest is maintained for between 20% to 50% of the mapping of the center of k-space. This, translates into correlating a period of substantially elevated or maximum rate of injection with the period of collection of image data corresponding to the center of k-space so that during at least 20% of the time of mapping k-space, a substantially elevated or maximum concentration of contrast agent is maintained in the artery of interest relative to adjacent veins; and preferably about 50%.

When imaging smaller arteries, in a preferred embodiment, an elevated or maximum concentration of contrast agent in the artery of interest is provided for greater than 50% of the time during which image data corresponding to the center of k-space is collected; and in a more preferred embodiment, an elevated or maximum concentration of contrast in the artery of interest is maintained for greater 75% of the mapping of the center of k-space. As a result, where the artery of interest is relatively small, the administration of the contrast agent may include a maximum or elevated rate of injection of the contrast agent of greater than 50% of the time of mapping of the center of k-space; and preferably between 50% to 85%, and most preferably greater than 75%. Under this circumstance, fewer artifacts are observed in the smaller vessels or arteries when the contrast is administered at a maximum or elevated rate over a longer period of the k-space mapping.

It Is noted that, in the aforementioned embodiment, the period of a maximum or substantially elevated rate of injection may also be compensated according to other factors including, the time delay due to the delivery system, the location of the artery, and the physical condition of the patient.

During the acquisition of image data corresponding to the center of k-space, it may be important to avoid excessively rapid changes in arterial contrast concentration. Rapidly changing blood signal during acquisition of the center of k-space may create image reconstruction artifacts. These artifacts may be minimized when the arterial signal intensity is uniform. Further, these artifacts may be minimized by slowly changing the arterial concentration during acquisition of image data and especially during acquisition of the center of k-space.

In those instances where the invention is implemented using paramagnetic contrast agents, infusion is at a rate that will provide a concentration of the agent in the arteries, such that the arteries will appear at least 50% brighter than any background structures, including veins, in the final image. In a preferred embodiment, the concentration of contrast agent will cause the longitudinal relaxation time (T1) of the protons in the arteries to be shorter than protons in any of the background material. Where the contrast agent causes the arteries to appear black in the final image (e.g., where the contrast agent shortens T2*, for example, Fe powder), the contrast agent should be infused at a rate and amount to insure that the effective transverse relaxation time (T2*) in the arteries is shorter than in any of the background material.

Any apparatus suitable for magnetic resonance imaging (MRI) of a portion of an animal body, for example, a human, may be used for acquisition of image data in the method of this invention. In particular, apparatus and imaging methods for magnetic resonance angiography are known in the art (see, e.g., U.S. Pat. Nos. 4,718,424; 5,034,694; and 5,167,232, incorporated herein by reference) and these may be used with the method of MRA with dynamic intravenous injection of magnetic resonance contrast agents taught herein, subject only to the constraints taught below.

The parameters of the imaging method of the present invention are discussed immediately below with respect to gadolinium chelates. It should be noted that other magnetic resonance contrast agents may be employed in practicing the present invention including paramagnetic contrast agents, such as those described by Marchal, et al., in Potchen, et al., eds., supra, pp. 305-322, the text of which is incorporated herein by reference.

Injection Parameters

Gadolinium chelates are paramagnetic agents which shorten the longitudinal relaxation time, T1, of blood according to EQUATION 1: ##EQU1## where: (1) the longitudinal relaxation time (T1) of blood without gadolinium is 1200 milliseconds; and

(2) Gd! is the blood concentration of a gadolinium chelate.

With reference to EQUATION 1, to achieve an arterial blood (T1) that is short compared to adjacent fat (T1=270) , it is necessary to substantially elevate the arterial blood concentration of the contrast agent in the artery of interest to be greater than (1/270 milliseconds-1/1200 milliseconds)/relaxivity of the contrast agent (or 2.9/seconds*relaxivity). Thus, the artery of interest includes a substantially elevated concentration of the contrast agent when that concentration is greater than 2.9 seconds⁻¹ relaxivity⁻¹ of the contrast agent.

A substantially elevated rate of infusion provides a substantially elevated concentration of the contrast agent in the artery of interest. That is, a substantially elevated rate of infusion provides an arterial blood concentration of the contrast in the artery of interest which is greater than 2.9 seconds⁻¹ relaxivity⁻¹ of the contrast.

As reflected in EQUATION 2, below, the arterial blood Gd! may be expressed in terms of the intravenous injection rate and the cardiac output during dynamic imaging at times short as compared to the recirculation time. ##EQU2##

As long as the gadolinium chelate is sufficiently small, the gadolinium chelate will rapidly redistribute into the extracellular compartment as it passes through the capillary bed and the venous concentration will be low or negligible compared to the arterial concentration. The relationship between the longitudinal relaxation time of arterial blood and the injection rate may then be determined by combining EQUATION 1 and EQUATION 2, as stated below in EQUATION 3: ##EQU3##

To achieve contrast between arterial blood and background tissue, the longitudinal relaxation time of the arterial blood should be reduced to less than that of the background tissues. Of all types of background tissues, fat (T1=270 msec) typically has the shortest longitudinal relaxation time. Assuming a typical minimum resting cardiac output of 0.0005 Liters/Kg-sec and requiring the longitudinal relaxation time to be less than 270 milliseconds simplifies EQUATION 3 to EQUATION 4 as shown below: ##EQU4##

By way of example, gadopentetate dimeglumine, gadodiamide, and gadoteridol are three paramagnetic gadolinium chelates that are readily available and rapidly redistribute into the extracellular fluid compartment. The relaxivities of gadopentetate dimeglumine and gadoteridol are 0.0045/molar-second. Based upon the aforementioned and using EQUATION 4, the minimum injection rate is greater than 0.033 millimole/Kg-minute.

With continued reference to EQUATION 4, a rate of infusion which is greater than 0.0015 Liters/Kg-sec² divided by the relaxivity may provide a maximum concentration of the contrast agent in the artery of interest. That is, infusing the contrast into the patient at a rate of greater than 0.0015 Liters/Kg-sec² divided by the relaxivity may yield a maximum arterial blood concentration of the paramagnetic contrast agent.

The total dose of gadolinium chelate required may be determined by multiplying the injection rate by the imaging time. For a relaxivity of 4.5/millimolar-second, and an imaging time of 5 minutes (300 seconds), the dose should substantially exceed 0.1 millimole/kilogram.

The dose of the gadolinium chelate may be within the range of 0.05 millimoles/kilogram body weight to 1 millimoles/kilogram body weight depending upon the time required to obtain the image. It should be noted that the dose of the contrast should not be too high such that there may be undesirable toxicity or T2 effects. In a preferred embodiment, the dose of the gadolinium chelate is within the range of 0.2 millimoles/kilogram body weight to 0.4 millimoles/kilogram body weight. In a more preferred embodiment, the dose of the gadolinium chelate is about 0.3 millimoles/kilogram body weight.

In those instances where the contrast injection times are longer than the recirculation time, the longitudinal relaxation time of arterial blood tends to be even shorter since a fraction of the gadolinium chelate will recirculate. It should be noted that a T1 of 270 ms (corresponding to the brightest background tissue fat) is equivalent to a gadopentetate dimeglumine concentration of 0.6 millimole/liter.

FIG. 1 illustrates the longitudinal relaxation time (T1) of blood as a function of infusion time and the total paramagnetic contrast dose for a paramagnetic contrast compound having a relativity of 4.5/millimolar-second. An examination of FIG. 1 reveals that the shortest T1 occurs with the shortest infusion time and the largest gadolinium dose. For typical imaging times of 3 to 5 minutes, FIG. 1 further reveals that the dose should be of the order of 0.2 millimoles/kilogram or larger in order to achieve a longitudinal relaxation time of blood significantly shorter than that of the brightest background tissue fat (T1=270) for the entire duration of imaging.

Imaging parameters

Any suitable T1 weighted magnetic resonance imaging sequence may be used during injection of the paramagnetic contrast. Suitable imaging sequences will be readily apparent to the skilled practitioner and are described in Potchen, et al., eds., supra. The following criteria for selection of preferred imaging parameters are based on experience in over 100 patients on a 1.5 Tesla General Electric signa magnet with version 4.7 software. A three-dimensional Fourier Transform (volume) acquisition (3D FT) is preferred in the abdomen because of its intrinsically high spatial resolution and high signal-to-noise ratio, even with a large, body coil. The gradient echo (gradient recalled) pulse sequences are preferred since they allow a short TR (repetition time) which allows a shorter imaging acquisition time. Short imaging times have the advantage of allowing the same total gadolinium dose to be injected at a faster rate.

Spoiled versus non-spoiled gradient echo imaging

It should be noted that one might expect steady state gradient echo imaging (GRASS) to be preferable to the spoiled gradient echo imaging because the long T2 (transverse relaxation time) of blood increases the steady state blood signal. However, this effect enhances veins more than arteries, because the fast, pulsatile flow of arterial blood spoils its steady state component. In theory, this may have the paradoxical effect of reduced arterial contrast. In practice, there may only be a small difference between the spoiled and unspoiled techniques. In patients with slow arterial flow (which is not self-spoiling), steady state gradient echo may be preferred. A spoiled gradient echo pulse sequence (SPGR) was chosen for most of the studies described herein to simplify the theory and analysis to reduce the potential for differential steady state magnetization between arterial blood, slower venous blood and background tissue.

Echo Time

Because the brightest background tissue is fat, it is preferable to use a TE (echo time) where fat and water are out of phase, thereby achieving an incremental improvement in vessel-to-background contrast. At 1.5 Tesla, this occurs about every 4.6 msec beginning at 2.3 msec which corresponds to a TE of 2.3, 6.9, 11.5, . . . msec. The shortest of these possible TE values (6.9 msec in the studies described herein) is preferred. Shorter TE's tend to minimize the effects of motion related phase dispersion.

Repetition Time

In a preferred embodiment, TR should be as short as is possible. A TR of 24-25 msec was the shortest possible on the equipment used for the studies described herein. As the TR is shortened, the flip angle must be adjusted to maintain the optimal T1 weighing.

Flip Angle

With a gadolinium chelate dose of 0.2 millimoles/kilogram and a 3-5 minute injection time and imaging time, the longitudinal relaxation time of the arterial blood is predicted to be in the order of 150 to 200 milliseconds. It will, however, be shorter as a result of the recirculation time being less than 3-5 minutes. The relative signal intensity, SI, in a 3D FT spoiled gradient echo acquisition as a function of blood T1, TR, T2, T2*, flip angle σ, and proton density N(H) may be expressed as stated in EQUATION 5, below, and calculated accordingly. ##EQU5##

FIG. 2 graphically illustrates relative signal intensity for T1 equal to 50, 100, 150, 270 (fat), and 1200 (blood) under the following conditions: (1) TR=25 milliseconds, and assuming TE is small compared to T2* (the observed transverse relaxation time) FIG. 2 reveals that a flip angle of about 40 degrees is optimal for maximizing blood-to-background tissue (fat) contrast when the longitudinal relaxation time (T1) of blood is of the order of 200 milliseconds. For larger gadolinium doses with faster injection rates, a larger flip angle may be more appropriate.

Volume Orientation

In order to minimize the image acquisition time, the imaging volume should be made as thin as possible while containing the arteries of interest. In this regard, it may be useful to orient the image volume for maximum in-plane coverage of the vessels of interest as opposed to the perpendicular orientation required for optimal time-of-flight magnetic resonance angiography. Optimizing the orientation and minimizing the thickness of the imaging volume is facilitated by first acquiring a conventional black-blood or time-of-flight MRI to use as a guide for precise localization. Phase and frequency encoding axes should be oriented such that cardiac and respiratory motion artifacts do not superimpose on the vessels of interest. Generally, for imaging the aorta-iliac system, the imaging volume should be oriented coronally, and the phase encoding axis should be set right-to-left. For imaging the thoracic aorta, a sagittal orientation is preferred and for imaging the subclavian arteries, an axial orientation is preferred.

Partitions

The number of partitions (slices) is determined by the thickness of the image volume divided by the partition thickness. The partition thickness is the image resolution along the axis perpendicular to the plane of the partitions. It may be useful to employ thin partitions in order to have high image resolution. The image acquisition time, however, linearly increases with the number of partitions. As a result, keeping the image acquisition time short requires minimizing the number of partitions.

It should be noted that there may be a loss of signal-to-noise as the voxel size is decreased by using higher resolution pixels. Generally, 0.5 to 2 millimeter resolution with 28 to 60 partitions is adequate for the aorta and major branch vessels. The skilled practitioner will balance the need to increase resolution by decreasing voxel size with the need to avoid excessive periods of time to acquire image data.

Field-of-View

The field-of-view must be made large enough to avoid excessive wrap-around artifact. Wrap around artifacts occur when there are structures outside the field of view along the phase encoding axis. These structures are mapped by the phase encoding process to superimpose on structures within the field of view.

In addition, because of the limited number of pixels along the frequency encoding axis and time penalty for each additional pixel along the phase encoding axis, it is also desirable to make the field-of-view as small as possible in order to maximize image resolution with the minimum image acquisition time. Generally, for imaging the abdominal or thoracic aorta, a field-of-view of about 36 centimeters is appropriate for most patients. It may be increased for larger patients and reduced for smaller patients. Smaller field-of-views may be used for other parts of the body.

Use of a no-phase wrap algorithm is a less preferred embodiment. Under the circumstance of this invention, this has a disadvantage of generally requiring more imaging time and, as a result, a larger gadolinium dose.

Coils

It is preferable to use the smallest possible coil in order to minimize noise. There is also an advantage to coils that encircle the body part of interest such that the signal will be homogeneous throughout the entire field-of-view.

Patient Positioning

The patient should be positioned such that the body part being imaged remains stationary during the acquisition of image.

Cardiac and Respiratory Motion Compensation

The phase artifact related to respiratory and cardiac motion may be minimized by combining the T1 weighted imaging sequence with respiratory or electrocardiographic gating. Gating has the disadvantage of increasing the scan time--particularly in patients with irregular rhythms. Compensation techniques in which the acquisition of the image data in k-space is matched to the respiratory and or cardiac cycle may eliminate some phase artifact without significantly increasing the scan time.

In imaging regions of the body that move substantially with respiration (e.g., the renal arteries) it may be useful to acquire data while the patient is holding his breath. This may require shortening the duration of the image acquisition time to under one minute. If the patient cannot hold his breath for the entire period of image acquisition, than it may be useful to hold the breath during acquisition of image data corresponding to the center of k-space and breathing only during acquisition of data corresponding to the periphery of the center of k-space.

Pre-scanning

The pre-scanning process is used to tune to the optimum frequency and to optimize the receiver gain. In the pre-scanning process, it is necessary to compensate for the changes in the patient's magnetic resonance signal that will occur during the contrast injection. In those instances when the paramagnetic contrast agent is a gadolinium chelate, it is preferable to tune to the water peak. About a 20% to 50% margin should be incorporated into the receiver gain setting to allow for increased signal during contrast administration.

Premedication

Premedicating patients with an analgesic or sedative such as diazepam may be useful for at least two reasons. Firstly, it may help the patient to tolerate the claustrophobic sensation of being within the magnet thereby reducing voluntary motion artifacts. Secondly, and more importantly, its relaxation and cardiac depressant effects tend to reduce the cardiac output. A lower cardiac output results in a higher arterial contrast concentration which thereby improves the image quality. This result is opposite from conventional magnetic resonance angiography which is degraded when the cardiac output decreases. By reducing the cardiac and respiratory rates, analgesics and sedatives may minimize the fraction of the image acquisition that is adversely affected by cardiac and respiratory motion artifacts.

Magnetic Resonance Contrast Agents

As mentioned above, many different magnetic resonance contrast agents may be employed when implementing the present invention; for example, numerous paramagnetic contrast agents are suitable. As mentioned above, gadolinium compounds, for example, paramagnetic gadolinium chelates, such as gadopentetate dimeglumine and gadoteridol, are readily available and rapidly redistribute into the extracellular fluid compartment. Other gadolinium compounds are acceptable. In general, preferred contrast agents have a high relaxivity, rapid redistribution into the extracellular fluid compartment, and are readily extracted from the capillary bed. It should be noted that, contrast agents that are extracted or degrade in the capillary bed are preferred in the present invention.

In particular, gadolinium chelates are commercially available from such companies as Bristol Meyers (under the name "Prohance") Berlex (under the name "Magnevist"), and Sterling Pharmaceutical (under the name "OmniScan"). It should be noted that the gadolinium chelate which is commercially available from Sterling Pharmaceutical, now Nycomed, appears to facilitate greater contrast enhancement between the artery and the surrounding veins and tissue.

Injection

In a preferred embodiment, the type or form of injection of the paramagnetic contrast is intravenous. The injection of the paramagnetic contrast is performed intravenously in order to eliminate or reduce the complications associated with the catheterization required for arterial injections.

The specific site of injection is important for several reasons. The site of injection should be remote from the "region of interest"; that is, the region that is to be scanned. For example, when imaging the abdominal aorta, intravenous injection of the paramagnetic contrast into an arm vein is preferred. Use of a leg vein should be avoided. Further, there may be some benefit in avoiding the antecubital fossa because the patient may bend the elbow during a long (3-5 minute) period of injection which may result in extravasation of the contrast into the subcutaneous tissues. As a result, under this condition, a forearm or upper arm vein may be preferable.

In those instances where an artery in the arm is to be imaged, the site of the injection may be a leg vein or a vein in the opposite arm. Here, the site of injection is remote from the "region of interest", i e., the artery in the arm.

Moreover, it is important to adapt the timing of a maximum or elevated rate of infusion of the contrast agent to correlate with the collection of image data which corresponds to the center of k-space. As will be discussed in more detail below, correlating the infusion rate with the mapping of k-space insures that the image data representative of the center of k-space is collected over some period during which a maximum or elevated concentration of contrast agent is maintained in the artery of interest relative to adjacent veins. There are several manners of adapting the timing of a maximum or substantially elevated rate of infusion to correlate a maximum or substantially elevated concentration of contrast agent with the mapping of k-space, including adjusting the timing to account for the time delay in the delivery or infusion apparatus, the size and location of the artery of interest, the type of pulse sequence employed by the imaging apparatus, and the condition of the patient.

In a preferred embodiment, as illustrated in FIGS. 3 and 4, a mechanical infusion or injection device 12 is an automated type of injector having reliable and consistent operating conditions. The infusion device 12 is employed to inject the magnetic resonance contrast agent into the vein of the patient at an infusion rate sufficient to provide contrast enhancement of an image of an artery relative to veins in the field of view of the magnetic resonance image and substantially throughout the period of acquisition of image data. The infusion device 12 couples to the patient using conventional techniques, for example, appropriately selected tubing 14 which permits fluid flow between the mechanical infusion device 12 and the patient. Such tubing may be, for example, an angiocatheter.

A mechanical injector is preferred because of the greater reliability and consistency when compared to injecting by hand. Since the magnetic field interferes with normal functioning of electronic devices, a pneumatic powered, spring loaded or other non-electric pump may be suitable. It should be noted, however, that an electrical pump may be used if its operation is unaffected by the operation of the magnetic resonance imaging system, e.g., if the pump is adequately shielded or if it is located sufficiently far from the magnet.

In one preferred embodiment, the mechanical infusion device 12 is coupled to the magnetic resonance imaging system 16 to facilitate proper or desired timing between the injection of the magnetic resonance contrast agent and the acquisition of the image data, in addition to providing proper or desired rates of infusion of the contrast agent.

In another preferred embodiment, proper or desired timing and rates of infusion of the contrast agent are controlled through a control mechanism at the mechanical infusion device 12. That is, the mechanism that controls the infusion timing and rates of infusion is implemented within the mechanical infusion device 12. In this circumstance, the mechanical infusion device 12 is a "self-contained" unit. For example, the infusion rate may be controlled with an adjustable fluid flow resistor.

As mentioned above, the infusion device 12 injects the magnetic resonance contrast in a controlled manner. The contrast may be contained in a vessel. As illustrated in FIGS. 3 and 4, the mechanical infusion device 12 is coupled to a vessel 18 which contains the magnetic resonance contrast agent. In one embodiment, the vessel 18 may contain a sufficient quantity of contrast agent for one application of the invention, e.g., a single use vessel. In an alternative embodiment, the vessel 18 may contain a quantity which allows several applications of the invention, e.g., a reservoir type of vessel. As is illustrated in FIG. 3, the mechanical infusion device 12 may be adapted to receive the vessel 18 somewhat like a fountain pen receiving an ink cartridge. In an alternative embodiment, as illustrated in FIG. 4, the infusion device 12 may be coupled to the vessel 18 using conventional methods.

FIGS. 5A and 5B illustrate a mechanical infusion device 12 in more detail. The mechanical infusion device 12 of FIGS. 5A and 5B include several of the components described and illustrated in U.S. Pat. Nos. 4,202,333; 4,298,000; 4,430,079; and 4,597,754. The descriptions of these patents, including several of the components of the mechanical infusion device 12 described therein, are incorporated herein by reference. Moreover, such infusion devices are commercially available from 3M Corporation and their product specification sheets are also incorporated by reference.

In those instances where the mechanical infusion device 12 is employed within the environment of the magnetic field, the infusion device 12 should be fitted or manufactured with magnetic resonance compatible material. For example, the infusion devices which are commercially available from 3M Corp., should be fitted with a magnetic resonance compatible spring. This requires manufacturing the spring from non-magnetic materials, for example, plastic or certain metal alloys such as eljaloy or inconel.

To implement the techniques of the present invention and thereby obtain a constant or variable rate of infusion of the magnetic resonance contrast agent, the device 12 of FIGS. 5A and 5B may include a spring which has a constant width and thickness in order to exert a constant force; or, alternatively, the spring may have a variable width and/or variable thickness to provide a variable spring force. Under this circumstance, the infusion rate may be controlled to be either constant or variable by design of the spring and, in effect, pre-programmed by selection of the spring's design parameters.

In one embodiment, the infusion device 12 infuses the magnetic resonance contrast agent at a variable infusion rate having a maximum or elevated rate of infusion which temporally correlates with the acquisition of image data which is representative of the center of k-space. Under this circumstance, data representative of the center of k-space is acquired during the period of a maximum or substantially elevated arterial gadolinium concentration in the artery of interest.

As mentioned above, in one preferred embodiment, the timing of a maximum or substantially elevated rate of injection may be correlated to the mapping of the center of k-space in accordance with the particulars of the delivery systems (e.g., physical length of the catheter/tubing), the location of the artery of interest, the size of the artery of interest, and/or the physical condition of the patient. Such a correlation technique allows image data corresponding to the center of k-space to be collected during a period of maximum or substantially elevated contrast concentration in the artery of interest.

It should be noted that the center of k-space may be characterized as 10% to 75% of the total k-space data which corresponds to the lowest spatial frequency information.

It should be further noted that a substantially elevated concentration of the contrast agent in the arterial blood may be described as a concentration which is greater than 2.9 seconds⁻¹ -relaxivity⁻¹. As mentioned above, a substantially elevated rate of infusion provides a substantially elevated concentration of the contrast agent in the artery of interest. In this regard, a substantially elevated rate of infusion provides a concentration of the contrast in the artery of interest which is greater than 2.9/seconds-relaxivity (of the contrast).

Further, as mentioned above, a rate of infusion of the contrast agent which is greater than 0.0015 Liters/Kg-sec² divided by the relaxivity may provide a maximum arterial blood concentration of the paramagnetic contrast agent in the artery of interest.

In one preferred embodiment, the infusion device 12 may be designed to accommodate a 50 cc syringe having a fluid capacity of 60 cc and containing one dose of the contrast agent. The infusion device 12 may also be designed to permit an external (manual) force on the syringe to modify or customize the rate of infusion of the contrast agent. This external force is separate from the force of the spring of the infusion device 12.

Further, the infusion device 12 may include a flow rate indicator 50 (FIG. 5B) to provide an indication of flow rate of the contrast agent to the patient. Under this circumstance, the operator may visually or audibly observe, in a rather simple manner, the rate of flow of the contrast agent. This will allow the operator to exert an external force more accurately (both in the amount of force applied and in a timing sense) thereby facilitating a modification of the predetermined injection rate. Further, in an automated-type of infusion device, the flow rate indicator 50 permits the operator to visually or audibly monitor a "pre-programmed" infusion rate or sequence.

Briefly, reference to FIG. 5A, the mechanical infusion device 12 further includes syringe 24, a syringe clamp 26, a syringe restraint or stop 28, a block and spring housing 30a, roller bearings 30b, a reflux valve 32 and a catheter 34 having tubing 34a and a needle 34b (butterfly type) . The syringe 24 contains the contrast agent to be administered to the patient during magnetic resonance imaging. A plunger 24a of the syringe 24 is engaged by a spring 40 which is housed in the block and spring housing 30a. In operation, the spring 40 engages the plunger 24a to pressurize the syringe 24. The syringe is maintained in a stationary position within the mechanical infusion device 12, and in particular, in housing base 12a, via the syringe clamp 26 and the syringe restraint 28.

In one preferred embodiment, the mechanical infusion device 12 is coupled to a saline drip apparatus 42 (saline drip 42a, tubing 42b and roller clamp 42c). The saline drip apparatus 42 is applied to an input of a y-port connector 44. The syringe 24 is applied to the other input of the y-port connector 44. This conventional configuration facilitates a saline flush following the administration of the contrast agent within the syringe 24. In those instances where the tubing leading from the saline drip apparatus 42 to the y-port connector 44 has a one-way valve to prevent reflux of contrast, it is acceptable to leave the saline drip "on" during infusion. Under this circumstance, as soon as the infusion of the contrast agent is complete, the drip infusion will automatically resume to "flush" gadolinium within the intravenous tubing and deliver the contrast agent which remains in the tubing to the patient.

The rate of injection of the contrast agent from the syringe 24 is determined or controlled, in large part, by the size or gauge of the needle 34b, which functions as a fluid flow restrictor according to Poiseulle's Law. The rate of injection is also controlled by the amount of force that the spring 40 (the restoring force of the spring 40) applies to the plunger 24a of syringe 24, the syringe cross-sectional area, the gadolinium viscosity as follows: ##EQU6## where: r=radius of flow restricting needle lumen;

F=spring force;

L=length of flow restricting needle;

μ=viscosity of the fluid; and

A=cross-sectional area of the syringe.

Examining the infusion rate equation immediately above reveals that a variation of the syringe size (A), needle length (L), and/or fluid viscosity (μ) impacts the rate of infusion of the contrast agent. The viscosity of the fluid, however, may be dependent on the temperature of the contrast agent (gadolinium chelate). Thus, in those instances where the temperature of the contrast agent alters the viscosity of the contrast, the rate of infusion is also dependent on this "variable."

It should be noted, however, that the influence of viscosity on the flow rate may be substantially reduced by employing a fluid flow restrictor which minimizes the effects of viscosity on the rate of fluid flow.

In a preferred embodiment, the characteristics of the spring 40 (e.g., spring force) may be selected or designed such that the spring 40 applies a constant force upon plunger 24a throughout the period of contrast infusion. In another preferred embodiment, the characteristics of the spring 40 may be selected or designed such that the spring 40 applies a variable force on the plunger 24a. That variable force may correlate with the imaging process so that a maximum or substantially elevated injection rate provides a maximum or substantially elevated concentration of contrast in the artery of interest during the collection of image data which corresponds to the center of k-space.

The rate of injection, however, may be increased or decreased using a manual, spring loaded, or pneumatic injection rate adjustment mechanism which may be connected to various components of the device 12, including, the spring 40, the block and spring housing 30a, the roller bearings 30b, the plunger 24a, the tubing 34a, and/or the fluid flow restrictor 34b. FIG. 5B illustrates a manual injection rate mechanism 50 for allowing the operator to readily alter the rate of injection and thereby modify the rate of injection of the contrast agent to accommodate or implement a desired timing of an elevated or maximum rate of flow of the contrast agent.

The spring force should be sufficient such that the flow restrictor, required to give the desired flow rate, has a flow resistance that is much greater than any flow resistance in the intravenous line. The spring force should not be so high that a person of ordinary strength can not reduce or increase this force when a manual spring adjustment mechanism 50 is designed as the means for adjusting the rate of flow (i.e., the amount of external force applied to the spring) of the contrast agent. In general, a spring with about 5-10 pounds of spring force is suitable.

A fluid flow restrictor may be manufactured from, include or be comprised of a needle, a short piece of tubing of narrow calibre (e.g. an intravenous angiocatheter of 20 gauge or larger may be satisfactory), an orifice (for example made of ruby or sapphire), a focal compression of the IV tubing, or other mechanism which impedes the flow of fluid.

It should be noted that a precision orifice may offer several advantages when employed as a fluid flow restrictor. For example, in those instances where an incompressible fluid is to be administered, such as gadopentetate dimeglumine, gadoteridol, or gadodiamide, flow through an orifice is governed by the Bernoulli effect. In this regard, the flow rate of the fluid through the orifice is proportional to the square root of the pressure drop: ##EQU7## where: K=a constant determined by the geometry of the orifice;

F=spring force; and

A=syringe cross-sectional area.

Further, it should be noted that the pressure drop across an orifice is governed by inertial effects of the fluid; the viscosity of the fluid has little to no impact. As a result, an orifice minimizes the influence of the viscosity of the fluid on the rate of flow of the fluid. Under this circumstance, by using an orifice as a fluid flow restrictor, the Bernoulli effect predicts the same flow rate regardless of temperature of the fluid and regardless of which gadolinium compound is employed. Although in practice it is essentially impossible to entirely eliminate viscosity effects of the fluid, those effect are markedly reduced.

TABLE 1 provides the infusion rate, with respect to three MR contrast agents, for a variety of needles and flow restricting orifices when employed in an infusion device 12 substantially as illustrated in FIG. 5A where the spring 40 is a 6 pound-force spring and the syringe 24 is a 1 inch diameter, 50 cc syringe.

In those instances where the rate of flow of the fluid is dependent on the ambient temperature or the temperature of the contrast agent, consistent operation of the infusion device 12 may require either a temperature controlled operating environment or use of a fluid flow restrictor whose operational characteristics are essentially unaffected by the viscosity of the fluid (e.g., a precision orifice).

With reference to FIG. 5B, in a preferred embodiment, the rate adjustment mechanism 50 is a manual type including a lever 50a by which the user may increase or decrease the force applied to the plunger 24a. The lever 50a engages the plunger 24a and spring 40 so that the resulting force applied to the plunger 24a is essentially determined by the sum of the force applied to the plunger 24a (i.e., by the lever 50a) and the spring force, F. By employing this configuration, the user may increase or decrease the rate of injection at a particular moment of the imaging sequence. For example, increasing the infusion rate at about 10 to about 40 seconds prior to the acquisition of image data corresponding to the center of the k-space would cause an elevated or relatively high arterial gadolinium level to be maintained in the artery of interest during acquisition of image data corresponding to the center of k-space (typically it takes about 10-40 seconds for venous blood in the arm to circulate through the heart and lungs to reach the artery of interest). Such a technique may provide additional contrast enhancement of the image of the artery relative to veins and surrounding tissue.

The infusion device 12 of FIG. 5B further includes a flow rate indicator 52 to provide the operator an indication of a flow rate (injection rate) of the contrast agent to the patient. Here, the operator may visually or audibly observe the rate of flow of the contrast agent to thereby accurately control the rate of injection of the contrast agent into the patient; the operator may customize or modify the contrast injection rate.

The flow rate indicator may be implemented using an optical type sensor for sensing the linear motion of, for example, the plunger 24a, the spring 40, and/or the block and spring housing 30a, or the rotational motion of the roller bearings 30b. Such a mechanism permits an accurate measurement with little to no impact on the operation of the injection device 12, including the motion of the plunger 24a and the operation of the spring 40. That is, an optical type rate indicator has an advantage of not requiring physical contact with the contrast agent in the syringe 24 or spring 40.

It is noted, however, that a fluid flow or motion sensor may also be employed in the flow rate indicator 50. Such devices provide accurate information regarding the rate of flow of the contrast agent in the syringe 24 or in the tubing 34a.

As mentioned above, when the mechanical infusion device 12 is employed within the environment of the magnetic field, the materials used to fabricate the device 12 should be non-magnetic. That is, magnetic materials should be avoided when the device 12 is implemented in or near the magnetic field of the magnetic resonance imaging apparatus. In those instances, the spring 40 (FIGS. 5A and 53) should be manufactured from non-magnetic materials, for example, eljaloy or inconel.

With reference to FIGS. 6A-C, the mechanical infusion device 12 may be implemented using a bag-cassette configuration. The bag 46 contains a contrast agent. Analogous to the syringe configuration of FIGS. 5A and 5E, the bag 46 may be placed into a cassette 48 which applies even pressure over the contact surface of the bag 46. In operation, the contrast agent then flows, similar to the syringe 24, from the bag, through the catheter 34 to the patient. As with the case with the syringe configuration, fluid flow control is provided by means of a fluid flow restrictor (i.e., the needle 34b) used in combination with a cassette 48 which provides the force.

It should be noted that the bag-cassette arrangement of FIGS. 6A-C may employ a saline drip apparatus 42 as well as a rate adjustment mechanism 50. As with the syringe configuration, the rate of injection may be increased or decreased using a manual, spring loaded, or pneumatic rate adjustment mechanism 50.

In some magnetic resonance suites, an opening exists in the wall dividing the magnet of the imaging apparatus and the control equipment (i.e., computer and other electronic devices). In these situations, standard infusion pumps (containing metal, magnetized material and electronic circuits) can be used from outside of the MR suite to implement the methods described herein.

In one preferred embodiment, a pump manufactured by Abbott, the Life Care 5000, may be implemented. The Life Care 5000 draws drugs (e.g., contrast agent) directly from a bottle and preloads it into a long length of tubing. The operating parameters of the Life Care 5000 may be preprogrammed to execute numerous infusion rates.

In another preferred embodiment, the injection rate for contrast is matched with the mapping of k-space so that a maximum or substantially elevated arterial gadolinium concentration correlates with acquisition of image data corresponding to the center of k-space. That is, the operating parameters of the pump may be pre-programmed to provide an injection rate for contrast agent which is matched with the mapping of k-space so that a maximum or substantially elevated rate of infusion occurs about 10-40 seconds prior to the collection of image data corresponding to the center of k-space.

This type of configuration offers several advantages including: (1) the contrast agent (gadolinium) need not be removed from its shipping containers into an intermediate container, for example, a syringe; (2) the programmability of the pump allows variable injection rates providing for a maximum rate at the peak when the center of k-space is being mapped (which may be the most critical period during image acquisition) ; (3) operator control of the operating parameters.

It should be noted that the Life Care 5000 Pump may not be ideally suited for implementing all of the techniques described herein. For example, such deficiencies include the rates of injection of the pump, the degree of programmability of the flow delivery characteristics of the pump, and allowing the pump to administer contrast from multiple containers which will permit multiple 20 cc vials to be used.

As mentioned above, the present invention controls the precise timing of the injection of the paramagnetic contrast to provide preferential enhancement of an artery of interest relative to adjacent veins and background tissue. The injection of the paramagnetic contrast agent should temporally correlate to a period during which imaging data is being collected and preferably concurrent with the acquisition of image data. When implementing longer pulse sequences (greater than 2 minutes) or pulse sequences which collect image data representative of the center of k-space some time after initiation of image data collection, it is important that no contrast be administered prior to magnetic resonance scan since the contrast may leak into the background tissues and cause a degradation of the image. If some paramagnetic contrast or other magnetic resonance contrast has been administered prior to imaging, it is preferred to delay the arterial scan until this contrast has been excreted by the patient, in order to increase the probability of obtaining the optimal images.

In those instances where the imaging apparatus employs pulse sequences having very short data acquisition periods the contrast agent may be injected before the initiation of collecting image data. Short pulse sequences may be characterized as those sequences for which the time required for contrast to circulate from injection site to the artery of interest becomes a significant fraction of the imaging time, for example, data acquisition periods of less than 2 minutes. Under this circumstance, injection of the contrast agent before acquisition of image data is necessary to allow circulation of the contrast agent in the patient and thereby correlate a maximum or substantially elevated arterial concentration with the collection of image data representing the center of k-space. Administering the contrast agent prior to the acquisition of image data would cause a relatively high arterial gadolinium level during the mapping of k-space. As discussed above, the relative timing between the administration of the contrast agent and the collection of image data representing the center of k-space should be adapted to account for the injection mechanism employed, the location of the artery of interest, the size of the artery of interest, and the physical condition of the patient. For example, the contrast may be administered about 10-40 seconds before collection of image data to account for venous blood in the arm to circulate through the heart and lungs to reach the artery of interest. Thus, the amount of time before acquisition of image data may depend on the configuration of the contrast delivery mechanism, the relative location of the artery of interest, the relative size of the artery of interest, and the condition of the patient, including the age of the patient. Employing these considerations in selecting and controlling the timing of the injection provide a more accurate alignment between the acquisition of data representative of the center of k-space and a period of maximum or substantially elevated contrast concentration in the artery of interest relative to adjacent veins.

In a preferred embodiment, a constant infusion should begin within a few seconds of initiation of the scan process. The contrast infusion should end about 20 or more seconds before the completion of the scan; this allows the intravenously injected contrast to circulate through the heart and into the arteries. A chaser of normal saline or other fluid may be used to insure injection of the entire dose of the paramagnetic contrast (e.g., gadolinium) and, in addition, to insure that there is sufficient venous return to propel the injected contrast to the heart. In a preferred embodiment, the contrast infusion rate is matched with the mapping of k-space so that the maximum arterial gadolinium concentration occurs during acquisition of the center of k-space. This may permit injecting over a shorter period of time to achieve either a higher injection rate or a lower contrast dose.

As discussed above, in order to adapt the timing of a maximum or substantially elevated rate of infusion and, in effect, a maximum or elevated arterial concentration of the contrast agent in the artery of interest, with the collection of image data corresponding to the k-space consideration should be given to the impact of the contrast agent delivery system or apparatus, the time interval for venous blood at the site of infusion to reach the artery of interest, the size of the artery of interest, and the location of the artery of interest has on this correlation. Each of these "variables" impact on the time between the maximum or elevated rate of injection and when a maximum or elevated arterial concentration in the artery of interest is observed.

The inherent delay in the contrast delivery system depends on, for example, the length of the tubing 14, between the infusion device 12 and the patient. Such a delay, once measured or determined, may be considered a "constant" for a given delivery or infusion apparatus.

The time interval for venous blood at the site of infusion to reach the artery of interest (i.e., the "circulation time") depends on the location of the artery of interest and the physical condition of the patient. Typically, this time interval may vary from about 10 to about 30 seconds. The time interval may be longer if the artery of interest is far from the heart, for example, in the foot, or if the site of infusion is a peripheral vein such as in the hand.

Further, young, healthy patients with normal cardiac output tend to have a faster flow and a short circulation time, between 10-20 seconds (typically about 15 seconds). Older patients and patients with low cardiac output have a slower circulation times, between 20-30 seconds or more. The circulation time will also be longer for an IV in the hand compared to one in the upper arm. Central lines which inject directly into the SVC or right atrium will have the shortest circulation times.

In another preferred embodiment, the period of a maximum or substantially elevated rate of infusion of the magnetic resonance contrast agent to the patient is adapted according to the size of the artery of interest to correlate with the period of the collection of image data corresponding to the center of k-space. In this regard, where the artery of interest is relatively large (e.g., the aorta), a period of a substantially elevated or maximum injection rate may overlap for a smaller fraction of the period of collecting image data representative of the center of k-space than where the artery is relatively small (e.g., renal). For example, when imaging larger arteries, the administration of the contrast agent may include a period of a substantially elevated or maximum rate of contrast which provides a substantially elevated or maximum arterial concentration for less than 50% of the period during which the system collects image data corresponding to the center of k-space; and preferably between 20% to 50%. Where the artery of interest is relatively small, it is preferable that a period of maximum or substantially elevated rate of injection provide a maximum or substantially elevated concentration of the contrast in the artery of interest for more than 50% of the period of mapping the center of k-space; and preferably between 50% to 85%.

With reference to FIG. 7, the infusion of the magnetic resonance contrast agent may be by way of manual means. In this embodiment, a syringe 20, having needle 22, is coupled to a vessel 14 containing the magnetic resonance contrast agent. The vessel 14 is coupled to the patient using conventional techniques, for example, appropriately selected tubing 14 which permits fluid flow between the vessel 18 and the patient, for example, an angiocatheter.

When implementing the present invention via manual injection, i.e., injecting the magnetic resonance contrast agent by hand, in a preferred embodiment, the infusion "path" includes a fluid flow restrictor which adds resistance to the flow of gadolinium during administration into the body. It should be noted that a fluid flow restrictor may be, for example, a standard injection needle or small calibre angiocatheter. In FIG. 7, the fluid flow restrictor may be the needle 22 of syringe 20 and/or the angiocatheter 14. Use of small needles, short pieces of tubing of narrow calibre, an orifice, and/or small calibre angiocatheters may alleviate errors of injecting the contrast too rapidly and, as a result, depleting or running-out of contrast before completion of the scan or improperly correlating a maximum or elevated rate of infusion with the mapping of k-space. In a preferred embodiment, the needle size may be 22 gauge or smaller diameter (higher than or equal to 22 gauge) depending upon the viscosity of the contrast agent.

It may be convenient to preload the contrast into a vessel or length of tubing with luer lock or other appropriate connectors at each end of the tubing. It is then possible to use a single saline filled syringe to inject the contrast followed by a saline chaser without having to switch syringes or pumps. Saline is a preferred fluid to use as a chaser since it can be made isotonic with blood and is compatible with most intravenous fluids and pharmaceuticals that may already be flowing through a patient's IV line.

In a preferred embodiment, the contrast is infused slowly at the beginning and fastest in the middle of the acquisition. This type of injection pattern, based upon the fact that the contrast does somewhat contribute to venous and background tissue enhancement, avoids excessive contrast early in the acquisition.

Post-Processing

Post-processing of the scan data may be used. Maximum intensity projection (MIP) collapse images are useful for rapidly examining the entire arterial circulation within the region of interest. It may be useful to reformat and selectively collapse the data through the specific arteries of interest. Additional contrast may be obtained by digitally subtracting a pre-gadolinium acquisition from the dynamic gadolinium acquisition. Volume rendering may also be useful and is possible with high contrast volume data sets.

Immediately below are examples of results obtained from use of preferred embodiments of the present invention. The parameters of the examples are detailed therein.

EXAMPLE 1

Contrast between peripheral arteries and veins in images obtained by imaging dynamically during the administration of gadopentetate dimeglumine was investigated in sixteen patients referred for aorta-iliac magnetic resonance arteriography. These included 9 males and 7 females with a mean age of 72 ranging from 67 to 83. The indications for the study included hypertension (6), abdominal aortic aneurysm (AAA, 6) claudication (4) and renal failure (9). Some patients had more than one indication.

Parameters

All imaging was performed on a 1.5 Tesla superconducting magnet (General Electric Medical Systems, Milwaukee, Wis.) using the body coil and version 4.7 software. A 3D FT, coronal, spoiled, gradient echo volume was acquired centered on the mid-abdomen. The imaging parameters included: 12 cm volume with 60 partitions, 2 mm partition thickness, TR of 25 msec, a TE of 6.9 msec, a flip angle of 40°, first order flow compensation, 36 centimeters field of view, 256 by 192 matrix. The imaging time was 5 minutes and 10 seconds. Frequency was set superior to inferior so that phase artifact from diaphragmatic and cardiac motion would not superimpose on the abdominal aorta and IVC. When possible, phase artifact noise was minimized by excluding the heart and lungs entirely from the field of view. No saturation pulses were employed. The volume data were reformatted through vessels of interest and also displayed as maximum intensity projections.

Gadopentetate Dimeglumine Injection

After pre-scanning, venous access was obtained via a 22 gauge angiocatheter. A dynamic acquisition was then performed during hand injection of gadopentetate dimeglumine (Berlex Laboratories, Cedar Knoll, N.J.), 0.2 millimoles/kilogram. The injection was initiated within 5 seconds of initiating the image acquisition. The injection rate was constant (within the limitations of a hand injection) and timed to last until 10-20 seconds before completion of the scan. The injection included a 5 cc normal saline chaser to ensure injection of the entire gadopentetate dimeglumine dose. As a result, the gadopentetate dimeglumine ended approximately 30-40 seconds before completion of the scan and the saline chaser ended about 10-20 seconds before completion of the scan. In order to compare to the conventional, non-dynamic, gadolinium-enhanced MRA, a second, identical acquisition was then acquired without altering the imaging or prescan parameters.

Signal Measurements

Signal intensity was measured in the abdominal aorta, IVC, iliac artery and vein, renal artery and vein, celiac trunk, SMA, portal vein, hepatic vein and background tissue (including fat, skeletal muscle, kidney, liver and spleen) for 7 regions of interest per measurement. As many of these measurements as possible were obtained from the central 20 partitions and all measurements were obtained from the central 40 partitions. Identical regions of interest were used to compare vessels on the dynamic and post-gadolinium images. The standard deviation of the aorta signal was recorded as noise. Differences in the aorta and IVC signal-to-noise ratio were evaluated for each patient as well as for the means of all patients with Students t-test. In addition, the significance of differences in the mean portal vein, hepatic vein, renal vein and iliac vein signal compared to the IVC were evaluated with Students t-test. The presence of aneurysms, occlusions and stenoses (>50%) was noted on the individual dynamic images and on maximum intensity projections and compared to findings at surgery or arteriography when available.

Results

All sixteen patients tolerated the imaging and gadopentetate dimeglumine well; there were no complications. FIGS. 8A-C illustrate the typical images obtained before, during and after injection of gadopentetate dimeglumine, respectively. Before the injection, the vessels were heavily saturated with only a few streaks of vessels visible at the edges of the 3D volume. Images obtained during injection showed enhancement of the arteries while the IVC remained indistinguishable from the background tissue. Aorta IVC signal intensity ratios, shown in TABLE 2, confirmed this preferential arterial enhancement in every patient studied. Images obtained after the injection was completed demonstrated comparable enhancement of both arteries and veins.

It should be noted that with dynamic imaging there is bright arterial as well as portal vein and splenic vein enhancement but no visible IVC or iliac vein enhancement and no in-plane saturation. Post gadopentetate dimeglumine images show comparable enhancement of both arteries and veins.

TABLE 3 provides the average signal intensity for all tissues studied for both the dynamic and post-injection sequences. With dynamic gadopentetate dimeglumine the average aorta signal-to-noise ratio was 10±0.9 compared to 5.1±1.4 in the IVC (p value≦0.0001), while post gadopentetate dimeglumine the aorta and IVC were nearly identical, 10±1.4 and 9.5±1.3 respectively. Although all veins were less bright than the aorta on the dynamic images compared to post gadopentetate dimeglumine images, there were variations among the veins analyzed. The iliac vein was the least enhanced, 4.7±1.6, while the portal vein was the brightest, 8.3±1.6 followed by the hepatic, 7.5±2.1, and renal, 6.2±1.8, veins; these differences were significant to the p<0.01 level compared to the mean IVC signal-to-noise ratio.

Angiographic and/or surgical correlation was available in 6 of the 16 patients. In the vascular segments for which definitive correlation was available, magnetic resonance arteriography correctly identified 2 occlusions (1 common iliac and 1 renal artery), 10 stenoses (4 renal artery, 2 iliac artery, 2 distal aorta, 1 inferior mesenteric artery and 1 celiac) and 6 aneurysms (3 aortic and 3 iliac artery). There was no evidence of arterial in-plane saturation in any patient. In one patient with a common iliac artery occlusion, there was no difficulty visualizing reconstituted flow distal to the occlusion.

TABLE 4 reveals an apparent trend for patients with a history of cardiac disease, claudication or aneurysms to have the greatest aorta/IVC signal intensity ratio. The sample size may have been too small to establish statistically significant correlations. Further, one patient with cardiac disease, aneurysmal disease and claudication had the highest aorta/IVC signal intensity ratio. These trends are opposite from time-of-flight imaging where cardiac disease and aneurysms are associated with image degradation.

EXAMPLE 2

In order to determine the optimal timing of contrast administration, two methods of dynamic administration, bolus and continuous infusion, were compared to non-dynamic injections and to conventional time-of-flight imaging.

Gadolinium enhanced magnetic resonance arteriography was performed in 52 patients referred for routine MRA of the abdominal aorta or branch vessels. Imaging was performed as described in Example 1. The total acquisition time was 5:08 minutes to cover approximately 36 cm of aorta and iliac artery in the superior to inferior dimension. In 20 of the patients of the dynamic gadolinium infusion imaging was performed with 28 partitions each 2 mm thick with a 256 by 256 matrix to reduce the scan time to 3:18 minutes.

After pre-scanning, venous access was obtained via a 22 gauge angiocatheter. A dynamic acquisition was then performed during hand injection of gadolinium dimeglumine (Berlex Laboratories, Cedar Knoll, N.J.) 0.2 millimoles/Kg. In 12 patients, the injection was given as a bolus. The bolus was begun within 5 seconds of starting the acquisition and completed within the first 1 to 2 minutes of the 5 minute scan. In the other 40 patients, an injection of the same dose was carefully timed to be constant and continuous over the entire period of imaging beginning within 5 seconds of commencing the acquisition and ending 20 seconds before the end of the acquisition. In all patients, a 5 cc normal saline chaser was given to ensure injection of the entire gadopentetate dimeglumine dose.

For comparison purposes, 16 of these patients were imaged with an identical acquisition after the dynamic infusion without altering the imaging or prescan parameters. Also, for comparison, axial 2D and multiple overlapping 3D (MOTSA) time-of-flight images were acquired prior to the gadolinium injection. Inferior presaturation pulses were used with the time-of-flight sequences to suppress venous in-flow.

Signal intensity was measured in all patients in the aorta, IVC and background tissues (fat and skeletal muscle) for at least 3 regions of interest per measurement for all sequences. The signal's standard deviation within the aorta was recorded as noise.

Images obtained dynamically, during steady infusion of gadopentetate dimeglumine, showed sufficient arterial enhancement to clearly define the aorta and branch vessel anatomy while the IVC and iliac veins remained indistinguishable from the background tissues. The portal vein is visible but is not as bright as the aorta. Images obtained non-dynamically, after the injection was completed, or with a dynamic bolus injection demonstrated comparable enhancements of both arteries and veins.

The observation of significant, preferential arterial enhancement with a continuous dynamic contrast infusion was confirmed by region of interest analysis (see TABLE 5 and FIG. 9). The ratio of aorta to IVC signal intensity for the 5 minute infusion, 2.0±0.5, was significantly higher than for non-dynamic imaging 1.1±0.1 (p<0.001) or for the dynamic bolus 1.2±0.2 (p<0.001). Even better differentiation between the aorta and IVC was obtained by injecting the same dose of gadopentetate dimeglumine more quickly over a 3:18 minute acquisition. Although this aorta-to-IVC signal intensity ratio was not as favorable as for 2D time-of-flight or MOTSA imaging, it was adequate in all cases for clearly distinguishing the aorta and abdominal aorta branch vessels from the IVC and iliac veins.

Dynamic contrast enhanced 3D imaging had no saturation, pulsatility or misregistration artifacts. Even in aneurysms, which tend to have stagnant and/or turbulent flow, there was no loss of signal. By comparison, every 2D time-of-flight study had pulsatility artifacts and some had misregistration and/or in-plane saturation artifacts. The MOTSA images had no pulsatility or misregistration artifacts but every MOTSA study showed some degree of arterial saturation and they were particularly degraded in aneurysmal segments.

Administering gadopentetate dimeglumine dynamically as a steady, continuous, infusion for the entire period of a 3D FT acquisition, at a dose of 0.2 millimoles/Kg, gives sufficient preferential arterial enhancement to visualize arteries distinctly from veins and background tissues regardless of the magnitude or direction of flow. The importance of injecting dynamically and continuously during the entire scan is illustrated by the absence of significant preferential enhancement when the contrast is administered non-dynamically or as a dynamic bolus. Images obtained at a lower dose, 0.1 millimole/Kg, were not useful.

Since dynamic gadolinium enhanced MRA does not depend upon the in-flow of unsaturated spins, it eliminates some of the saturation problems that complicate routine time-of-flight imaging. The imaging volume can be oriented in any plane for optimal coverage of the vessels of interest without concern for saturation. In these patients, in-plane, coronal imaging of the aorta-iliac system reduced the image acquisition time by 5 to 20 fold over 2D time-of-flight and MOTSA imaging and had superior resolution and superior aorta signal-to-noise ratios.

A 3D FT acquisition was used in this example partly because of its intrinsically high spatial resolution and high signal-to-noise and also because its main limitation, arterial saturation, is eliminated by the gadolinium. The TE was chosen to be as short as possible at a value where fat and water protons are out of phase. A short TE helps to minimize motion related phase artifacts. Having fat and water out of phase provides an element of fat suppression which improves artery-to-background contrast since fat is the brightest background tissue.

EXAMPLE 3

MRA image data for a patient presenting with an abdominal aortic aneurysm was acquired as described in Example 1. MRA images are shown in FIGS. 10A and 10B.

The MRA of FIG. 10A depicts the aneurysmal aorta and aneurysmal common iliac arteries as well as severe stenoses of the right external iliac (curved arrow) and inferior mesenteric (straight arrow) arteries and a mild stenosis of the left common iliac artery. The internal iliac arteries are excluded because of their posterior course. FIG. 10B illustrates a digital subtraction angiogram which confirms the findings in FIG. 10A as discussed immediately above.

EXAMPLE 4

A pump (as illustrated in FIG. 5A) was loaded with a 50 cc syringe containing 42 cc of gadodiamide. A 23 gauge butterfly was attached to the end of the syringe with its standard luer-lock connector and plugged into a side port of the patient's intravenous (IV) line within a few feet of the IV skin entry site. The pump was located approximately 15 cm or more away from the imaged volume.

Usually the IV site was in the forearm or antecubital fossa. A plastic Kelly clamp on the butterfly tubing prevented premature gadolinium infusion. The gadolinium infusion was begun simultaneous with beginning the image acquisition by releasing the clamp on the butterfly tubing. This combination of 6 pounds spring force and a 5/8 inch, 23 gauge butterfly needle gave an infusion rate of 18 cc/minute which was slightly reduced by the additional resistance of the IV tubing and angiocatheter. For the 42 cc volume of gadolinium, the calculated infusion time was 2:20 minutes. This was shortened by manually increasing the rate of injection during the middle of the acquisition such that the maximum arterial concentration occurred during acquisition of the center of k-space.

The pump infusion finished with one minute of scan time remaining. Residual gadolinium within the IV tubing (about 4 cc) was flushed through with saline to ensure delivery or the entire dose.

The procedure of EXAMPLE 4 produced excellent quality MRA images of arteries without the confounding effects of excessive venous enhancement.

Various preferred embodiments of the present invention have been described. It is understood, however, that changes, modifications and permutations can be made without departing from the true scope and spirit of the present invention as defined by the following claims, which are to be interpreted in view of the foregoing.

For example, the correlation techniques for adapting the timing of a maximum or elevated rate of infusion and the collection of image data representative of the center of k-space may be used in various combinations according to the particulars of the contrast delivery system, the artery of interest, the location of the artery of interest, the size of the artery of interest, the condition of the patient, and the type of pulse sequence employed by the imaging apparatus. Adapting the timing of a maximum or elevated rate of infusion to correlate with the collection of image data representative of the center of k-space insures that such data will be collected over some period during which a maximum or elevated contrast concentration is realized in the artery of interest.

                  TABLE 1                                                          ______________________________________                                         Infusion Rates of Gadolinium Chelates at 24° C.                                      Gd-DTPA Gadoteridol                                                                              Gadodiamide                                     ______________________________________                                         Viscosity @ 20° C.  cP!*                                                               4.9       2.0       2.0                                         ______________________________________                                         Flow Restrictor                                                                               Infusion Rate                                                   Needles  Size      Gd-DTPA   Gadoteridol                                                                            Gadodiamide                               ______________________________________                                         BD ® 18g 1.5"  100       126     120                                       Terumo ®                                                                            20g 1.5"  44        66      64                                        Terumo ®                                                                            21g 1.5"  27        47      44                                        Terumo ®                                                                            22g 1.5"  15        29      23                                        Terumo ®                                                                            23g 1"    <4        <4      <4                                        ______________________________________                                         Butterflies                                                                             Size      Gd-DTPA   Gadoteridol                                                                            Gadodiamide                               ______________________________________                                         ABBOTT ®                                                                            21g .75"  21        37      36                                        ABBOTT ®                                                                            23g .75"   8        19      18                                        ABBOTT ®                                                                            25g .375" <4         8      7.3                                       ______________________________________                                         Orifice  Size      Gd-DTPA   Gadoteridol                                                                            Gadodiamide                               ______________________________________                                                  0.010"    21        25      25                                        ______________________________________                                          *Values provided by manufacturer (Nycomed)                               

                                      TABLE 2                                      __________________________________________________________________________     Aorta/IVC Signal Intensity Ratios for Dynamic 3D Imaging                       Patient                                                                               Primary                                                                              heart    Signal Intensity During Injection                        #-sex                                                                              Age                                                                               Indication                                                                           disease                                                                            creatinine                                                                          Aorta                                                                               IVC  ratio**                                                                            p value                                    __________________________________________________________________________     1-m 83 AAA   yes*                                                                               2    7.9 ± 1.0                                                                        3.9 ± 0.6                                                                        2.0 <.0001                                     2-f 73 hypertension                                                                         yes*                                                                               .8   11 ± 1.0                                                                         8.2 ± 1.3                                                                        1.4 .0002                                      3-m 73 claudication                                                                         yes 2.2  10 ± 2.0                                                                         3.7 ± 0.5                                                                        2.8 .0003                                      4-f 67 hypertension                                                                         no  .9   10 ± 0.4                                                                         5.1 ± 0.6                                                                        2.0 <.0001                                     5-f 70 hypertension                                                                         yes*                                                                               3    8.9 ± 0.9                                                                        4.5 ± 0.4                                                                        2.0 <.0001                                     6-m 67 renal failure                                                                        yes 6    11 ± 0.5                                                                         4.9 ± 0.4                                                                        2.2 <.0001                                     7-f 80 AAA   yes 1.8  10 ± 0.4                                                                         5.9 ± 0.5                                                                        1.8 <.0001                                     8-f 76 renal failure                                                                        yes*                                                                               3.6  9.1 ± 0.6                                                                        5.0 ± 0.6                                                                        1.8 <.0001                                     9-m 68 AAA   no  1    11 ± 0.5                                                                         7.2 ± 0.3                                                                        1.4 <.0001                                     10-m                                                                               70 claudication                                                                         yes 1.2  11 ± 0.5                                                                         5.4 ± 0.3                                                                        2.0 <.0001                                     11-m                                                                               74 hypertension                                                                         no  1    8.9 ± 0.3                                                                        6.0 ± 0.8                                                                        1.5 <.0001                                     12-m                                                                               80 hypertension                                                                         yes*                                                                               3.2  10 ± 0.4                                                                         3.8 ± 0.9                                                                        2.6 <.0001                                     13-m                                                                               74 AAA   yes 4    9.8 ± 1.0                                                                        3.7 ± 0.8                                                                        2.6 <.0001                                     14-f                                                                               67 AAA   no  1    10 ± 0.3                                                                         5.9 ± 0.6                                                                        1.8 <.0001                                     15-m                                                                               67 hypertension                                                                         no  1.5  11 ± 0.9                                                                         4.6 ± 0.9                                                                        2.4 <.0001                                     16-f                                                                               71 claudication                                                                         yes*                                                                               6    11 ± 1.3                                                                         3.1 ± 0.6                                                                        3.5 <.0001                                     AVERAGE               10 ± 0.9                                                                         5.1 ± 1.4                                                                        2.0 <0.0001                                    __________________________________________________________________________      *cardiac disease with history of CHF                                           **Aorta/IVC signal intensity ratio                                       

                  TABLE 3                                                          ______________________________________                                         Average Signal-To-Noise Ratios                                                 During and Post Gadopentetate Dimeglumine Injection                                      Dynamic Post      Ratio                                                        Injection                                                                              Injection Dynamic/Post                                       ______________________________________                                         ARTERIES                                                                       Aorta       10 ± 0.9                                                                              10 ± 1.4                                                                              1.0                                            Iliac Artery                                                                               9.8 ± 1.3                                                                             10 ± 1.3                                                                              .98                                            Renal Artery                                                                               9.7 ± 1.9                                                                             10 ± 2.5                                                                              .99                                            Celiac & SMA                                                                               10 ± 1.7                                                                              11 ± 1.8                                                                              .91                                            VEINS                                                                          IVC         5.1 ± 1.4                                                                             9.5 ± 1.3**                                                                           .54                                            Iliac Vein  4.7 ± 1.6*                                                                            9.2 ± 1.3**                                                                           .51                                            Renal Vein  6.2 ± 1.8*                                                                            9.1 ± 1.9**                                                                           .68                                            Hepatic Vein                                                                               7.5 ± 2.1*                                                                            8.3 ± 1.0**                                                                           .90                                            Portal Vein 8.3 ± 1.6*                                                                            9.0 ± 3.3**                                                                           .92                                            BACKGROUND                                                                     Kidney      7.3 ± 1.0                                                                             8.3 ± 1.0                                                                             .88                                            Liver       5.3 ± 0.6                                                                             5.8 ± 1.8                                                                             .91                                            Spleen      5.9 ± 2.3                                                                             6.3 ± 2.3                                                                             1.1                                            Fat         4.3 ± 0.7                                                                             4.0 ± 0.8                                                                             1.1                                            Muscle      2.4 ± 0.5                                                                             3.2 ± 0.7                                                                             .75                                            ______________________________________                                          *p > 0.01 compared to IVC signal intensity                                     **p > 0.01 compared to signal intensity for dynamic injection                  ***standard deviation of signal in the space outside the patient         

                  TABLE 4                                                          ______________________________________                                         Effect of Cardiac Disease, Claudication and                                    Aneurysms on Aorta/IVC Signal Intensity Ratio                                  Subgroup     # of Patients                                                                              Aorta/IVC*                                                                               p value                                     ______________________________________                                         Cardiac Disease                                                                             12          2.2 ± 0.6                                                                             0.08                                        No Cardiac Disease                                                                          4           1.8 ± 0.4                                                                             0.08                                        Claudication 4           2.6 ± 0.8                                                                             0.12                                        No Claudication                                                                             12          2.0 ± 0.4                                                                             0.12                                        Aneurysm     7           2.2 ± 0.7                                                                             0.32                                        No Aneurysm  9           2.0 ± 0.5                                                                             0.32                                        ______________________________________                                          *Signal Intensity Ratio                                                  

                                      TABLE 5                                      __________________________________________________________________________     Effect of Injection Method on                                                  Aorta Signal-to-Noise and Contrast-to-Noise Ratios                                   Contrast   Image                                                                               Voxel                    Aorta-                                                                              Aorta-                                                                              Aorta-                Pulse Injection                                                                            # of time/cm                                                                             Volume                                                                              Saturation                                                                          Pulsatility                                                                         Aorta                                                                               Aorta/IVC                                                                           IVC  fat  muscle                Sequence                                                                             Method                                                                               patients                                                                            (sec/cm)                                                                            (mm.sup.3)                                                                          Artifacts                                                                           Artifacts                                                                           SNR  SI ratio                                                                            CNR  CNR  CNR                   __________________________________________________________________________     2D TOF                                                                               No gado                                                                              11   40   6.0  yes  yes  8.2 ± 2.8                                                                        3.7 ± 1                                                                          5.8 ± 1.9                                                                        5.5                                                                                 6.8 ± 2.4          MOTSA No gado                                                                              12   92   4.7  yes  no   8.9 ± 2.5                                                                        2.7 ± 0.9                                                                        5.1 ± 1.8                                                                        4.9                                                                                 6.3 ± 1.9          Gado: 3D                                                                             non-dynamic                                                                          16   9    3.1  no   no   9 ± 2.0                                                                          1.1 ± 0.1                                                                        0.6 ± 0.5                                                                        5.4                                                                                 6.2 ± 1.9          Gado: 3D                                                                             bolus*                                                                               12   9    3.1  no   no   12 ± 2.4                                                                         1.2 ± 0.2                                                                        2.7 ± 1.4                                                                        7.5                                                                                 9.1 ± 1.9          Gado: 3D                                                                             infusion**.sup.1                                                                     20   9    3.1  no   no   10 ± 1.2                                                                         2.0 ± 0.5                                                                        4.7 ± 1.4                                                                        5.4                                                                                 7.3 ± 1.1          Gado: 3D                                                                             infusion**.sup.2                                                                     20   5.5  3.1  no   no   10 ± 2                                                                           2.4 ± 0.8                                                                        5.6 ± 1.7                                                                        6.8                                                                                 8.2 ±              __________________________________________________________________________                                                              1.7                    SNR = signalto-noise ratio                                                     CNR = contrastto-noise ratio                                                   *gadopentetate dimeglumine given dynamically as a bolus within the first       minutes of the acquisition.                                                    **gadopentetate dimeglumine given dynamically as a constant infusion           spread over the entire acquisition.                                            .sup.1 5 minutes                                                               .sup.2 3 minutes                                                          

What is claimed is:
 1. A method of imaging a portion of an artery in a patient using magnetic resonance imaging, comprising,collecting a first image data set; administering a magnetic resonance contrast agent to the patient; and collecting a second image data set including collecting image data which is representative of a substantial portion of a center of k-space while a concentration of the contrast agent in the artery is substantially higher than a concentration of contrast agent in adjacent veins and background tissue.
 2. The method of claim 1 wherein the steps of collecting the first and second image data sets include using 3D gradient echo pulse sequences and the step of administering the magnetic resonance contrast agent includes administering the contrast agent to the patient by bolus injection.
 3. The method of claim 2 wherein at least one of the 3D gradient echo pulse sequences is a 3D spoiled gradient echo pulse sequence.
 4. The method of claim 1 wherein the step of collecting first image data set includes collecting data before administering a magnetic resonance contrast agent to the patient.
 5. The method of claim 1 wherein the step of collecting the second image data set further includes collecting image data which is representative of the center of k-space while a flip angle is between about 20° to about 90°, TR≦25 milliseconds and TE≦7 milliseconds.
 6. The method of claim 1 further including constructing an image of the artery by subtracting at least two images selected from a set of images comprising the images produced by the first image data set and images produced by the second image data set.
 7. The method of claim 6 wherein the step of collecting the second image data set further includes collecting at least 25% of the image data which is representative of the center of k-space while the concentration of contrast agent in the artery is substantially a maximum.
 8. The method of claim 6 wherein the step of collecting the second image data set further includes collecting image data which is representative of a substantial portion of the center of k-space while the concentration of contrast agent in the artery is substantially a maximum.
 9. The method of claim 6 wherein the step of administering the magnetic resonance contrast agent to the patient further includes administering the contrast agent to the patient substantially throughout collecting the second image data set.
 10. The method of claim 1 wherein the step of collecting the second image data set further includes collecting a portion of the image data which is representative of a center of k-space while the concentration of the contrast agent in the artery is substantially a maximum.
 11. A method of imaging a portion of an artery in a patient using magnetic resonance imaging, the method comprising,collecting a first image data set before administering a substantial amount of a magnetic resonance contrast agent; administering a magnetic resonance contrast agent to the patient; and collecting a second image data set including collecting image data which is representative of at least a portion of a center of k-space while a concentration of the contrast agent in the artery is about a maximum relative to a concentration of contrast agent in adjacent veins and background tissue.
 12. The method of claim 11 wherein the step of administering a magnetic resonance contrast agent to the patient includes administering the contrast agent to the patient by bolus injection.
 13. The method of claim 11 wherein the step of collecting the second image data set further includes collecting a substantial portion of the image data which is representative of the center of k-space while the concentration of contrast agent in the artery is substantially higher than a concentration of contrast agent in adjacent veins and background tissue.
 14. The method of claim 11 wherein the steps of collecting the first and second image data sets further include using 3D gradient echo pulse sequences.
 15. The method of claim 14 wherein at least one of 3D gradient echo pulse sequences is a 3D spoiled gradient echo pulse sequence.
 16. The method of claim 11 wherein the step of collecting the second image data set further includes collecting a substantial portion of the center of k-space while a flip angle is between about 20° to about 90°, TR≦30 milliseconds and TE≦7 milliseconds.
 17. The method of claim 11 wherein the step of collecting the second image data set further includes collecting the image data which is representative of the center of k-space during administration of the contrast agent and while the concentration of contrast agent in the artery is substantially higher than a concentration of contrast agent in adjacent veins and background tissue.
 18. The method of claim 11 wherein the step of collecting the second image data set further includes collecting substantially all of the center of k-space while the concentration of the contrast agent in the artery is substantially greater than a concentration of the contrast agent in veins adjacent to the artery.
 19. The method of claim 11 further includes constructing an image of the artery using the second image data set.
 20. The method of claim 19 further includes constructing an image of the artery by subtracting at least two images selected from a set of images comprising the images produced by the first image data set and images produced by the second image data set.
 21. A method of imaging a portion of an artery in a patient using magnetic resonance imaging, the portion of the artery being within a region of interest, the method comprising,collecting a first image data set administering a magnetic resonance contrast agent to the patient; and collecting a second image data set including:collecting image data which is representative of at least a portion of a center of k-space while a concentration of the contrast agent in the artery is about a maximum relative to a concentration of contrast agent in adjacent veins and background tissue; and collecting image data which is representative of a substantial portion of the center of k-space while a concentration of the contrast agent in the artery is substantially higher than a concentration of the contrast agent in veins and background tissue adjacent to the artery.
 22. The method of claim 21 wherein the step of administering a magnetic resonance contrast agent to the patient includes administering the contrast agent to the patient by bolus injection and the step of collecting the second image data set includes using a 3D gradient echo pulse sequence.
 23. The method of claim 22 wherein at least one of the 3D gradient echo pulse sequences is a 3D spoiled gradient echo pulse sequence.
 24. The method of claim 21 further includes constructing an image of the artery by subtracting at least two images selected from a set of images comprising the images produced by the first image data set and images produced by the second image data set.
 25. The method of claim 21 wherein the step of collecting the second image data set further includes collecting a substantial portion of the image data which is representative of the center of k-space while the concentration of contrast agent in the artery is substantially higher than a concentration of contrast agent in adjacent veins and background tissue and while a flip angle is between about 20° to about 90°, TR≦30 milliseconds and TE≦7 milliseconds.
 26. The method of claim 21 wherein the step of collecting image data which is representative of the substantial portion of the center of k-space includes collecting image data while a slice thickness is less than 5 millimeters.
 27. The method of claim 21 wherein the step of collecting image data which is representative of the substantial portion of the center of k-space includes collecting image data using a coronally oriented imaging volume.
 28. The method of claim 21 further including constructing an image of the artery by subtracting the image data of the first image data set from the image data of the second image data set.
 29. The method of claim 21 further including constructing an image of the artery by subtracting the image data of the second image data set from the image data of the first image data set.
 30. The method of claim 21 wherein the step of collecting image data which is representative of the substantial portion of the center of k-space includes collecting image data using a coronally oriented imaging volume and wherein a field of view is greater than about 30 centimeters.
 31. A method of imaging an artery in a patient using magnetic resonance imaging, comprising,collecting a first image data set; administering a magnetic resonance contrast agent to the patient; and collecting a second image data set using a body coil, the step of collecting the second image data set includes collecting image data which is representative of a substantial portion of a center of k-space while a concentration of the contrast agent in the artery is substantially higher than a concentration of the contrast agent in adjacent veins and background tissue.
 32. The method of claim 31 wherein the steps of collecting the first and second image data sets includes using 3D gradient echo pulse sequences.
 33. The method of claim 31 wherein at least one of the 3D gradient echo pulse sequences is a 3D spoiled gradient echo pulse sequence.
 34. The method of claim 31 wherein the step of collecting the second image data set further includes collecting a portion of the center of k-space while the concentration of the contrast agent in the artery is substantially a maximum.
 35. The method of claim 31 wherein the step of collecting the second image data set further includes collecting the image data which is representative of the center of k-space while a flip angle is between about 20° to about 90°, TR≦30 milliseconds and TE≦7 milliseconds.
 36. The method of claim 31 wherein the step of collecting the second image data set further includes collecting at least 50% of the image data which is representative of the center of k-space while the concentration of the contrast agent in the artery is substantially higher than the concentration of the contrast agent in veins adjacent to the artery and while a flip angle is between about 20° to about 90°, TR≦30 milliseconds and TE≦7 milliseconds.
 37. The method of claim 31 wherein the step of collecting the second image data set further includes collecting image data which is representative of the substantial portion of the center of k-space during administration of the contrast agent and while a flip angle is between about 25° to about 90°.
 38. The method of claim 31 further includes constructing an image of the artery using the second image data set.
 39. The method of claim 31 further includes constructing an image of the artery by subtracting at least two images selected from a set of images comprising the images produced by the first image data set and images produced by the second image data set.
 40. The method of claim 39 wherein the step of collecting the second image data set further includes collecting the image data which is representative of the center of k-space while a flip angle is between about 20° to about 90°, TR≦30 milliseconds and TE≦7 milliseconds.
 41. The method of claim 31 wherein the step of collecting image data which is representative of the center of k-space includes collecting image data while a slice thickness is less than 5 millimeters.
 42. The method of claim 31 wherein the step of collecting image data which is representative of the center of k-space includes collecting image data using a coronally oriented imaging volume and wherein a field of view is greater than about 30 centimeters.
 43. The method of claim 31 further including constructing an image of the artery by subtracting the image data of the first image data set from the image data of the second image data set.
 44. The method of claim 31 further including constructing an image of the artery by subtracting the image data of the second image data set from the image data of the first image data set. 